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Abstract 

It is commonly believed as a fundamental principle that energy-momentum conservation of a 
physical system is the result of space-time symmetry. However, for classical particle-field systems, 
e.g., charged particles interacting through self-consistent electromagnetic or electrostatic fields, 
such a connection has only been cautiously suggested. It has not been formally established. The 
difficulty is due to the fact that the dynamics of particles and the electromagnetic helds reside 
on different manifolds. We show how to overcome this difficulty and establish the connection 
by generalizing the Euler-Lagrange equation, the central component of a held theory, to a so- 
called weak form. The weak Euler-Lagrange equation induces a new type of hux, called the weak 
Euler-Lagrange current, which enters conservation laws. Using held theory together with the weak 
Euler-Lagrange equation developed here, energy-momentum conservation laws that are difficult to 
hnd otherwise can be systematically derived from the underlying space-time symmetry. 
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I. INTRODUCTION 


It has been widely accepted as a fundamental principle of physics that energy-momentum 
conservation of a classical or quantum system is due to the underlying space-time symmetry 
that the system admits. However, for classical systems with particles and self-generated 
interacting fields, the connection between ener gy- momentum conservation and space-time 
symmetry has only been cautiously suggested l| and has not been formally established. 
Examples of such classical particle-field systems include charged particles in an accelerator 
or a magnetic conhnement device interacting through the self-consistent electromagnetic 
helds j^. To understand and overcome this difficulty, we need to examine the details of 
the held theory. In the standard held theory, one writes down a Lagrangian density L, and 
the associated Euler-Lagrange equation determines the dynamics of the system. When the 
Euler-Lagrange equation is satished, a symmetry condition is equivalent to a conservation 
law. This is of course the celebrated Noether’s theorem 3|, 1^. It is surprising to hnd out 
for classical particle-held systems that the standard Euler-Lagrange equation does not hold 
anymore. This is because the dynamics of the particles and the electromagnetic helds reside 
on diherent manifolds. The electromagnetic helds are dehned on the space-time domain, 
whereas the particle trajectories as a held are only dehned on the time-axis. This is why the 
link between the symmetry and conservation law breaks down for these systems. This unique 
feature has not been discussed before, and it makes a signihcant diherence in the formulation 
of the held theory presented here. What we have discovered is that when the standard Euler- 
Lagrange equation breaks down, the held equations of these systems assume a more general 
form that can be viewed as a weak Euler-Lagrange equation. It a pleasant surprise to hnd 
out that this weak Euler-Lagrange equation can also link symmetries with conservation 
laws as in the standard held theory, where the regular Euler-Lagrange equation provides 
the link. The diherence is that the weak Euler-Lagrange equation induces a new type of 
current (unknown previously), called the weak Euler-Lagrange current, in conservation laws, 
in addition to the Noether current for the standard held theory. For many classical particle- 


,3 


or attracting 


held systems, such as particles interacting through electrostatic potentials 
Newtonian potentials 3 2|, energy-momentum conservation laws are difficult to hnd. Using 
the held theory with the weak Euler-Lagrange equation developed here, energy-momentum 
conservation laws can be systematically derived from the underlying space-time symmetries. 
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This paper is organized as follows. In Sec. m the classical particle-field systems and 
the difficnlty of establishing the connections between symmetries and conservation laws are 
introdnced. The weak Enler-Lagrange eqnation and its role in estabilishing conservation 
laws are given in Sec. lIIIl The last section snmmerizes the main resnlts of the paper. 


II. CLASSICAL PARTICLE-FIELD SYSTEMS 


The classical non-relativistic particle-field system in flat space is governed by the Newton- 
Maxwell eqnations 
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where Xsp{t) as a fnnction of time is the trajectory of the p-th particle of the s-species, 
and qs and rUg are the particle charge and mass, respectively. The electric field E{x,t) and 
the magnetic held B{x,t) are fnnctions of space-time. Eqnations ([I])-([2l can be expressed 
eqnivalently in the form of the Klimontovich-Maxwell (KM) eqnations 
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where Fs{x,v,t) = J2pS{Xsp — x)6{Xsp — v) is the Klimontovich distribntion fnnction in 
the phase space {x^v). 

Rednced models are often nsed in plasma physics. For example, the electrostatic 
Klimontovich-Poisson (KP) system is given by 
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where Bq{x) is a background magnetic field produced by steady external currents, and 
E = -V 6 is the longitudinal electric field. Another well-known reduced model is the 

fir 

Klimontovich-Darwin (KD) system |8l-lll|. 
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The local energy-momentum conservation laws for the Klimontovich-Maxwell system 


is well-known 
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where we have introduced 82 = 8 {Xgp — x) to simplify the notation. Through the following 
identities, 
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the conservation laws can be expressed equivalently in terms of the distribution function 
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For the reduced systems, e.g., the KP system and the KD system, it is also critical 
to know the exact local energy-momentum conservation laws admitted by the models. In 
practical applications, such as current drive and heating with lower-hybrid waves 1^, and 
electrostatic drift-wave turbulence, such local energy-momentum conservation laws for the 


reduced system have profound implications 
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the exact conservation laws admitted by the KP and KD systems, which are not exact special 
cases of the KM system, and shonld be viewed as independent systems in their own right. 
For example, we cannot take the exact energy-momentnm eqnations fllSp and (ITB]) . and 
approximate E by — and B by Bq to obtain the exact energy-momentnm conservation 
law for the KP system, even thongh the conservation law obtained this way conld be an 
approximate one for the KP system. The existence of exact local conservation laws is a 
necessary condition for the models to be theoretically well-posed and for the validity of 
particle simulations based on the KP or KD systems jlO[ |. 


On the other hand, conservation laws and symmetries are closely related. It is commonly 
believed that, according to Noether’s theorem 3|, 1^, conservation laws can be derived from 
the symmetries of the corresponding held theories. In standard held theories, this is cer¬ 
tainly true, and the symmetry in time for the action is related to energy conservation, and 
the symmetry in space corresponds to momentum conservation. Therefore, it is reasonable 
to expect that by analyzing the symmetries of the actions and Lagrangian densities for the 
reduced systems considered here, we may be able to systematically derive the desired con¬ 
servation laws. However, it is surprising to hnd out that for particle-held systems considered 
here, the held theory works diherently. First, let’s recall the action and Lagrangian density 
for the KM system given by Low H. 
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from 6 A/6Xsp = 0, 6 A/6(j) = 

0, and 


6 A/6A = 0. For the KP system, the action and Lagrangian density are given by 
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where Aq is the vector potential for a given external magnetic held Bq = V x Aq. For the 
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KD system, the action and Lagrangian density for the KM system are 
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Based on the spirit of Noether’s theorem, we would like to determine whether the local 
energy-momentum conservation laws can be derived from the symmetries of the correspond¬ 
ing Lagrangian density. It turns out that the answer to this question is not as simple as 
that in standard held theory. This is because the helds in the present held theory, i.e., 
Xsp{t), (j){x,t), and A{x,t) are dehned on diherent domains. The potentials are dehned 
on the space-time domain {x,t), whereas the particle trajectory Xsp{t) is only dehned on 
the time-axis. This unique feature has not been discussed before, and it makes a signihcant 
diherence in the formulation of the held theory presented here. 

In the next section, we develop the held theory for classical particle-systems with this fea¬ 
ture, in particular, the KM system, the KP system, and the KD system. The most distinct 
characteristic of the held theory presented here is that the held equation for Xsp{t) assumes 
a form we call the weak Euler-Lagrange (EL) equation, which is diherent from the standard 
Euler-Lagrange equation. The necessity of using the weak EL equation is mandated by the 
fact that Xsp{t), as a held, is not dehned on the entire space-time domain, but only on the 
time-axis. The weak EL equation with respect to Xsp{t) plays an indispensable role in the 
symmetry analysis and derivation of local conservation laws. For the KP system and KD sys¬ 
tem, the analysis developed here enables us to determine the desired local conservation laws, 
which have not been systematically discussed in the literature. For the KM system, where 
the local energy-momentum conservation laws ffT^ - flTBD are known, the present analysis 
serves the purpose of establishing a connection between the energy-momentum conservation 
laws and symmetries of the LaCTangian density Lkm- Interestingly, such a connection has 
only been cautiously suggested [ij but not explicitly established previously. This is perhaps 
not surprising, because the weak EL equation developed here is needed to establish the 
connection. Due to the space limitation, we present in Sec. lIIIl the detailed derivation of the 
held theory and new conservation laws only for the KP system of a magnetized plasma, and 
summarize the main results for the KM and the KD systems at the end. 
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In plasma physics, one often works with the Vlasov-Maxwell (VM) system. Equations 
dH])-® recover the VM equations when two-particle correlations (collisions) become negligi¬ 
bly small as the number of particles becomes increasingly large, while holding total charge 
and and charge to mass ratio hxed. In the present study, we work with the Klimontovich- 
Maxwell system (ED-® or ([I])-® and pass to the limit of the Vlasov-Maxwell system when 
necessary under the assumption of negligible collisions. Similarly, the Vlasov-Poisson (VP) 
and Vlasov-Darwin (VD) systems are regarded as the collisionless limits of the KP and KD 
systems, respectively. As a reduced system, the KP (or VP) system describes many impor¬ 
tant physical processes when the characteristic velocity of the particles or waves are much 
slower than the speed of light. These include electrostatic waves in plasmas (Langmuir 


0 . 


waves) 
moving with the beam 


and collective dynamics and excitations in charged particle beams in a frame 
The fundamental theory of Landau damping 17|] was hrst devel¬ 


oped for the VP (or KP) system. In astrophysics, the VP (or KP) system has also been used 


to model 
potential 


associated Landau d 
physics community 


;he collective dynamics of self-gravitating systems with an attractive Newtonian 
a, l7|. Because of these important applications, the VP (or KP) system and its 
am ping have also been studied with great interest in the mathematical 


18 


19j. We also note that while our focus here is on particle-held sys¬ 


tems, 
et al. 


erian held theories for the VM and VP systems have been developed by Morrison 
2[H23l| using a variety of theoretical constructions. In Eulerian theories, the particle 


distribution in phase space replaces Xsp{t) as the held variable. 


III. WEAK EULER-LAGRANGE EQUATION, SYMMETRY, AND CONSERVA¬ 
TION LAWS 


We begin with Eq. fl2T1) for the KP system, and determine how the action and Lagrangian 
density vary in response to the held variation SXgp and 6 (j){x,t), 
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In Eq. fl2^ . (p^i = dp/dx^ and integration by parts has been applied with respect to terms 
containing dExp/dp^i and dLxp/dXgp. Here, E^{Lkp) and Ex^piLxp) are the Euler op- 
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erators with respect to 0 and X^p, respectively. For a variable h, Dh/Dx"^ and Dh/Dt 
represent the space-time derivatives when h = h{x,t) is considered as a held on the space- 
time domain. Because 6 (j){x,t) is arbitrary, 6 A/6(j) = 0 requires the Euler-Lagrange (EL) 
equation for 0 to hold, i.e., E^{Lkp) = 0, which is indeed the Poisson equation flTH . as 
expected. The held equation for X^p is more interesting. Because SXgp is arbitrary only on 
the time-axis, the condition 5A/5Xsp = 0 requires only that the integral of Ex^piLxp) over 
space vanish, i.e.. 


(I^xExALkp) = 0. 


(31) 


Equation flTID will be called the submanifold Euler-Lagrangian equation because it is dehned 
only on the time-axis after the integrating over the spatial variable. If Xgp were a function 
of the entire space-time domain, then Ex^piLxp) would vanish everywhere, as in the case 
for 0(a;,f). In general, we expect that Ex^piL^p) 7^ 0. 

We now derive an explicit expression for Ex^ {L^p)- For the hrst term in Ex^ (L^p), 
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where the momentum Pgp density and Hamiltonian H^p density are dehned as 
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The second term in Ex^piLxp) is given by 
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Substituting Eq. (l35|) into the submanifold EL equation (l31l) . we immediately recover New¬ 
ton’s equation for Xgp, i.e., 
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which reduces Eq. fl55]) to 


Ex.,{Lkp) - (Hsp Xsp ■ Psp) + • (XspPsp) . (37) 

As expected, Ex^piLxp) ^ 0 . We will refer to Eq. flTTl) as the weak Euler-Lagrange equation, 
which is the foundation for the subsequent analysis of the local conservation laws. The 
qualiher “weak” is used to indicate the fact that only the spatial integral of the Euler 
derivative Ex^piLxp) is zero [see Eq. dTID ]. in comparison with the standard EL equation, 
which demands that the Euler derivative vanishes everywhere. 

We dehne a symmetry of the action A[(j), X^p] to be a group of transformation 

(s, t, 0, Xgpj I y (s?, f, 0, Xgpj (38) 


such that 
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If the symmetry is generated by a vector held on the space of {x, t, 0, Xgp), 
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then the inhnitesimal criteria of invariance is given by 




prV{L) + L Div^ = 0, 
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where Div^ is the divergence of the vector held 




(41) 


on the space-time domain, and prV is the prolongation of the vector held V on [x, t, 0, Xgp). 
The prolongation prV is a vector held on the jet space, consisting of the space of {x, t, 0, Xgp) 
and the space of derivatives of (0, J^) with respect to {x,t)- A comprehensive description 
of this subject can be found in Ref. |^. Given the symmetry vector held V, the inhnitesimal 
criteria for invariance will generate the desired conservation law corresponding to the sym¬ 
metry vector held V, after use is made of the EL equation as well as the weak EL equation 
for the systems in the present study. We hrst look for the symmetry group that generates 
local energy conservation. The group of transformation 


[x,t,(f>,Xgp) = {x,t + e,4),Xgp), eeR 


(42) 
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is a symmetry of Lkp, because Lkp does not depend on t explicitly, i.e., dExp/dt = 0, 
which can be written as 
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Equation (H3|l is the special form of Eq. (1401) for this symmetry group. From the EL equation 
for ((), i.e., E^{Lkp) = 0, we obtain 
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The weak EL equation for Xgp, i.e., Eq. flHTl) . gives 
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Combining Eqs. (|44|) and (|45|) . we obtain the first local energy conservation law. 


E 

m 


(v^ 

Stt 


s,p 


m 


^ ■V'2 



d 

j ^2 

dx 


-1 

dvr 


s,p 




We subtract the identify 
1 d r 


1 d 


4. at + 4^8^ ■ “ 

from Eq. fl46p to express the energy conservation law in another equivalent form 
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In terms of the distribution function Fg, we obtain 
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We emphasize again that Eq. 0491) is the exact energy conservation law admitted by the 
KP system Eqs. (jH]) and flTU]) . and it cannot be obtained by replacing E by — and B by 
Bo in the energy conservation law for the KM sytem (1151) . The sum of the last two terms 
in Eq. 0491) is the electrostatic Poynting flux of the KP system, first discussed by Similon 


24l | for an unmagnetized plasma by algebraic manipulation. Its importance for electrostatic 
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particle simulations was addressed by Decyk [25|. Here, it appears naturally as a conse¬ 
quence of the symmetry analysis. We observe that the external Bo does not contribute to 
the energy flux of the electromagnetic held. To further appreciate the importance of the 
exact energy conservation law flT^ . let’s consider the well-established technique of current 
drive and heating of a magnetized plasma using electrostatic lower-hybrid (LH) waves 12| , 
which are adequately described by the KM (or KP) system. In this application, the energy 
and momentum of the LH waves are converted to those of the particles, and it is of prac¬ 
tical importance to know the heating power of a specihc LH wave system. However, if we 
calculated the energy hux of the LH waves from Eq. flT^ by replacing E by — V0 and B by 
Bo, we would hnd that V • [V0 x Sq] = 0, i.e., the LH waves do not carry an energy hux. 
This is obviously erroneous. The typical power of such systems in modern magnetic fusion 
devices is several megawatts. The correct way to calculate the energy hux of the LH waves 
is to use Eq. flT^ instead. Specihcally, the several megawatts of energy carried by the LH 
waves how into the plasma through the last two terms in Eq. (09]). 

Up to now, we have treated the KM and KP systems as independent systems, each of 
which has its own governing equations, Lagrangian, and conservation laws. On the other 
hand, it is also correct to treat the KP system as the electrostatic approximation to the KM 
system. From the perspective of the governing equations, this approximation is equivalent 
to replacing E by —and B by Bo in the KM system. But this simple procedure does 
not work for the corresponding conservation laws. What is needed here is a more rigorous 
procedure to derive the electrostatic approximation that reduces from the KM system to 
the KP system. After this rigorous procedure is carried out, we hnd that the correct energy 
conservation law for the KP system obtained from that of the KM system is actually Eq. (091) . 
instead of that obtained from Eq. (1151) by replacing E by — V0 and B by Bo- A similar 
argument applies to the momentum conservation law for the KP system, i.e., Eq. (1571) . These 
derivations are given in detail in the Appendix. 

Our next goal is to search for the symmetry that generates the momentum conservation 
law. In standard held theories, if the Lagrangian density does not depend on x explicitly, 
then it admits the symmetry of spatial translation, x = x + eu, for a constant vector u and 
e E R. Then the usual form of Noether’s theorem leads to momentum conservation. This 
strategy does not work here because Lkp depends on x explicitly through 62 = S{Xsp — x) 
and Ao{x). However, if we simultaneously translate both x and X^p by the same amount. 
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then 62 is invariant. Thus, we consider the translational transformation 
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under which (j){x) = (j){x) = (j)[x — eu) and Xsp{i) = Xsp{t) + eu. When Ao{x) = 0, we can 
verify that Eq. (11(1 is satished, and Eq. flKl is indeed a symmtry admitted by Lkp- The 
corresponding vector field is 
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, and V is the only non-vanishing component of PrV since it is a constant. The notation 
d/dx here represents d/dx^ for i = 1,2, 3. In this case, the inhnitesimal criteria of invariance 
in (lid is 
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It will be clear shortly that this term represents a part of the momentum input due to the 
external magnetic held through the Lorentz force. For the hrst term in Eq. (I53|) . we invoke 
the EL equation E^{Lkp) = 0 to obtain 
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For the second term in Eq. (I5d . the weak EL equation for Xgp (llTIl is applied, which gives 
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Therefore, the conservation law generated by Eq. (I5ip is 
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Evidently, this is the local momentum conservation. In terms of the distribution function 
Fg, it can be expressed as 
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The first term on the left-hand side of Eq. 0561) or Eq. 0571) is the rate of variation of the 
momentum density, the second term is the divergence of the flux, and the term on the 
right-hand side is the momentum input due to the background magnetic held. Note that 
the momentum density is purely mechanical, and does not include the electromagnetic mo¬ 
mentum density — V0 x Bq/Attc. This is not totally intuitive. This conservation law is the 
result of the symmetry 05UI) . which is different from the well-known translational symmetry 
for standard held theory. Because Lkp depends on x explicitly through 62 = S{Xsp — x), 
a translation in x alone is not a symmetry of Lkp, even when Af){x) = 0. Instead, the 
symmetry group (1501) simultaneously translates the space x and the held Xgp by the same 
amount. 

For the KD system, the weak EL equation for Xgp is 
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Energy conservation follows from the inhnitesimal criteria (HOI) for the symmetry transfor¬ 
mation (H2|) after the weak EL equation (j58i) for Xgp and the EL equations for 0 and A are 
applied, i.e.. 
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Stt 


S ^ 




rrtsV 3 


d'^v 




rrisV 3 

Fs - vd V -I- 


At + E X B 


An 


= 0 . 


Similarly, the inhnitesimal criteria for the symmetry group (l5(J|) gives the momentum con¬ 
servation relation 

E X B\ 


d 

dt VV 


Fgvd'^v + 


47r 


+ 


A 

dx 


E 


m.o 


FgVvd^v + 


(V0)^ + B^ + 2V0 -EE + BB-AtAt 


8n 


An 


= 0 . 
(59) 


For the KM system, the weak EL equation for X^p is 
dLxD D dLxD 


^XspiLKo) = 


sp 


Dt dX 


sp 


E 

dx 


-A ■ Xsp + (j) - 2 ^sp ^ 


\ 1 

d 

r * / ■ \ 1 

)52 

dx 

[x,p + a) 82 


(60) 


The symmetry groups (142 p and (I5UI) gives the energy and momentum conservation laws ([T5p 
and (ITBp after the weak EL equation for Xgp and EL equations for 0 and A are applied. 
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IV. SUMMARY AND CONCLUSIONS 


In summary, a close examination of the field theory for classical particle-field systems re¬ 
veals that the particle field Xgp and the electromagnetic field reside on different manifolds. 
This unique feature is fount to imply that Ex^piL), the Euler derivative of the Lagrangian 
density L with respect to particle’s trajectory Xgp, does not vanish on the space-time man¬ 
ifold, which is surprisingly different from the standard field theory. In fact. 


,T\_ dL D ( dL \ 



( 61 ) 


for some non-vanishing tensor T. Equation fIM]) is what we call the weak Euler-Lagrange 
equation, and it is the most essential component in establishing the connection between 
energy-momentum conservation and space-time symmetry for classical particle-field systems. 
In fact, the energy-momentum conservation law follows from the infinitesimal criteria of the 
space-time system, after the weak Euler-Lagrange equation is applied. The non-vanishing 
tensor T is a new type of flux called the weak Euler-Lagrange current that enters the con¬ 
servation laws. For the Klimontovich-Maxwell (or Vlasov-Maxwell) system, this theoretical 
construction explicitly links the well-known energy-momentum conservation law with the 
space-time symmetry, which was only cautiously suggested previously. For reduced systems, 
such as the Klimontovich-Poisson (or Vlasov-Poisson) system and the Klimontovich-Darwin 
(Vlasov-Darwin) system, this theoretical construction enable us to start from fundamental 
symmetry properties in order to systematically derive the energy-momentum conservation 
laws, which are difficult to determine otherwise. 


APPENDIX 

In the Sec. m we have treated the Klimontovich-Maxwell (KM) and Klimontovich- 
Poisson (KP) systems as independent systems, each of which has its own governing equations, 
Lagrangian, and conservation laws. The Vlasov-Poisson (VP) and Vlasov-Darwin (VD) sys¬ 
tems are regarded as the collisionless limits of the KP and KD systems, respectively. As a 
reduced system, the KP (or VP) system describes many important physical processes when 
the characteristic velocity of the particles or waves are much slower than the speed of light. 
These include electrostatic waves in plasmas (Langmuir waves) j^, and collective dynamics 
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and excitations in charged particle beams in a frame moving with the beam 


mental theory of Landan damping 


S- 


The fnnda- 


17( 1 was hrst developed for the VP (or KP) system. In 


astrophysics, the VP (or KP) system has also been nsed to model the collective dynamics 


of self-gravitating systems with an attractive Newtonian potential 7|. Becanse of these 
important applications, the VP (or KP) system and its associated Landan dam ping have 
also been stndied with great interest in the mathematical physics commnnity 18l. Il9|. 


On the other hand, it is also correct to treat the KP system as the electrostatic approxima¬ 
tion to the KM system. From the perspective of the governing eqnations, this approximation 
is eqnivalent to replacing E by —V0 and B by Bq in the KM system. Bnt this simple pro- 
cednre does not work for the corresponding conservation laws. In this section, we present a 
more rigorons procednre to carry ont the electrostatic approximation that passes from the 
KM system to the KP system. After this rigorons procednre is carried ont, we hnd that the 
correct energy conservation law for the KP system obtained from that of the KM system is 
actually Eq. flT^ . instead of that obtained from Eq. f[T5]) by replacing E by — V0 and B by 
Bo. 


The electrostatic approximation applies when the characteristic velocity of the particles 
V and phase velocity of the waves u/k is much slower than the speed of light, i.e., when 
v/c^u/ckr.^e<^ 1. When this condition is satished, it turns out that the KM (or VP) 
system admits solutions with the following ordering 


El =Ef + eE\^'^ + e^E\^^ + 

(62) 

Et=^E'f^ + 0{e^), 

(63) 

B =Bo + + O(e^), 

(64) 

F, + e'Ff) + 0 {e^), 

(65) 


where Ei and Et are the longitudinal and transverse components of the electric held, respec¬ 
tively, and Bo is the externally applied magnetic held with V x Bo = 0 inside the plasma. 
The superscripts “(0)”, “(1)”, and “(2)’’represent the orders e°, e^, and e^. To the leading 
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order in e, i.e., 0(e°), the KM system is 


dFP 

dt 

(-V0 + YxBo) 
ox \mj s \ c J 

5ir(o) 

•^=0’ 

Ov 

(66) 


VV = -47rE?.y 


(67) 


O 

III 

1 

<1 


(68) 


which is indeed the KP system. Higher-order equations can be derived in a straightforward 
manner. 

For present purposes, we only need the hrst-order equation for the hrst-order magnetic 
held 


V X = — E / vF^dv + 
^ S J 


IdE 


( 0 ) 


dt 


(69) 


Note that is determined by the leading-order helds and due to the fact that 
v/c ^ oj/ck ~ e -C 1. Even though B^^^ does not enter the governing equations for the KP 
system, i.e., Eqs. (1661) and fl671) . it will enter the leading order energy conservation law for 
the KP system. Starting from the exact energy conservation law for the KM system, i.e., 
Eq. (8), we retain all the leading-order terms to obtain 




(v^ 

Svr 


+ E / 


m.v 


d^v 


+ V 


-cV0 X 
Air 


+ E / Ff’O^vd^v 


= 0. (70) 


The term —cV0 x B^^'^/Air is the Poynting hux due to the leading-order longitudinal electric 
held eI^^ and the hrst-order magnetic held B^^\ and it must be included in the leading- 
order energy equation, because ck/oj ~ 1/e raises the order of this term by one. Since B^^'> 
is uniquely determined by and Ef^'^ through Eq. fl69|) . the Poynting hux term can be 
expressed as 


V- 


-cV0 X B^F 


Air 


= — V0 • V X 
Air 


Y^qsj vF/^Uv - ^ . V0 = V Ye 5.0 y F/^^vd^v 


Air 


(t)V4>, 


t ) 


(71) 


where the continuity equation derived from Eq. 066 p has been used. Finally, the leading-order 
energy conservation equation is 


E 

m 

E_ 

dx 


(V0)^ 

Sir 


moV ,q 

-^—d^v 


+ 11 i 

S J 

E f F/^^'^vd^v + ^q^cj) J F/^^vd^v- 


Air 


0V0/ = 0, 


(72) 


16 
























which is identical to Eq. fll5|) . if is identified with Fg. This demonstrates that the energy 
conservation derived from the held theoretical approach for the KP system is not only more 
rigorous in mathematical treatment, but also more correct in physics content than the simple 
approach of replacing E by — V0 and B by Bq in the energy conservation law for the KM 
system. A similar argument applies to the momentum conservation law for the KP system, 
i.e., Eq. (lT6l) . 
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